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This work investigates the utility of femtosecond pulse shaping in increasing the efficiency of 
Raman excitation of molecules in the strong-field interaction regime. We study experimentally 
and theoretically the elTect of pulse shaping on the strength of non-resonant coherent anti-Stokes 
Raman scattering in iodine vapor at laser intensities exceeding lO'^^ W/cm^. We show that unlike the 
perturbative case, shaping strong non-resonant laser pulses can increase the signal strength beyond 
that observed with the transform-limited excitation. Both adiabatic and non-adiabatic schemes of 
excitation are explored, and the dilTerences of their potential in increasing the excitation efficiency 
are discussed. 
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INTRODUCTION 

The broad spectral bandwidth of ultrashort laser 
pulses is often used to excite coherent molecular wave 
packets consisting of a number of rotational and vibra- 
tional eigenstates. Studying the dynamics of an excited 
wave packet on femtosecond time scale represents a pop- 
ular approach to molecular spectroscopy [1]. Femtosec- 
ond spectroscopy benefits from the technology of pulse 
shaping [5] which offers selectivity and control of molec- 
ular excitation and has been successfully implemented 
in various applications such as multiphoton ionization 
[3] , stimulated Raman scattering [3l H] , four- wave mixing 
[3 [B] , coherent anti-Stokes Raman scattering [7HT^ and 
coherent control of chemical reactions [T^ . 

In the weak-field regime of Raman excitation of 
molecules, i.e. when molecular states are not changed 
significantly by the applied laser fields and the perturba- 
tion theory holds, the complex amplitudes of the excited 
molecular states are proportional to the corresponding 
resonant spectral components of the two-photon excita- 
tion field [m [16] . In the absence of intermediate res- 
onances, an upper bound on the absolute efficiency of 
exciting a particular mode is set by the available laser in- 
tensity and is reached with unshaped transform-limited 
(TL) laser pulses. Hence, improving the efficiency of an 
off-resonance Raman process in the perturbative regime 
can be achieved by increasing the laser intensity, whereas 
shaping the spectrum of the driving field can only sup- 
press rather than enhance the excitation at any given fre- 
quency. Stronger laser fields modify the molecular field- 
free states, most importantly via AC Stark shifts, often 
suppressing the rate of the target transition [TT'-'SO]. As 
a result, unshaped pulses no longer provide maximum ef- 
ficiency of transfering molecules to the target vibrational 
state in the strong- field limit. 

In this work, we investigate the utility of femtosecond 



pulse shaping to increase the magnitude of non-resonant 
vibrational excitation. We explore, both experimentally 
and theoretically, the efficiency of exciting vibrational 
wave packets in the ground electronic state of molecular 
iodine subject to strong laser pulses (> 10^'^ W/cm^). 
Evaluating the efficiency of strong-field induced vibra- 
tional excitation by means of coherent anti-Stokes Raman 
scattering (CARS), we demonstrate that pulse shaping 
can lead to the enhancement of the nonlinear spectro- 
scopic signal by at least 50%. 

Two qualitatively different approaches are analyzed. 
First, we suppress strong-field effects (such as AC Stark 
shifts) that may reduce the efficiency of the CARS pro- 
cess. In contrast to the feedback-loop based adaptive 
control [3 [T71 [2T], we achieve this by applying pre- 
determined pulse shapes to the excitation pulses in such 
a way as to lower the instantaneous field strength while 
preserving the amplitude of the two-photon field at the 
frequency of Raman resonances. 

In the second approach we try to employ the technique 
of adiabatic passage (AP) to enhance off-resonance vibra- 
tional excitation by making use of, rather than avoiding, 
strong-field effects. Non-adiabatic population transfer is 
very sensitive to the fluctuations of the laser parameters 
like intensity or pulse duration |22j . In AP, one exploits 
an adiabatic time evolution to improve the robustness 
of the process and enhance the efficiency of population 
transfer. Earlier works have demonstrated a number of 
successful implementations of adiabatic transfer in two- 
and three-level systems driven by resonant laser radi- 
ation [23H27] . Previous studies questioned the utility 
of AP in multi-level molecular systems interacting with 
strong femtosecond pulses [2H1 HH] ■ In this work we put 
to test the techniques of two-photon rapid adiabatic pas- 
sage (RAP (301 EJ), stimulated rapid adiabatic passage 
(STIRAP [23) and "piecewise AP" [32l|33]. The latter 
approach has been proven successful in atomic systems 
[271 [51] but has not been implemented with molecules. 



We use femtosecond CARS to determine the efRciency 
of vibrational excitation in the ground electronic state 
of molecular iodine. Two strong laser pulses, pump and 
Stokes, with frequencies Wp and cos ^-nd intensities ex- 
ceeding 10^'^ W/cm^, prepare the vibrational coherence 
via an off-resonance two-photon Raman transition. A 
weak probe pulse with frequency ujpr, separated in time 
from pump and Stokes, scatters off this coherence gen- 
erating the anti-Stokes field at tUaS — ~ + ^pr- 
The detected CARS signal at w^s serves as a quantita- 
tive measure of the degree of vibrational excitations, and 
its dependence on the spectral shapes of both pump and 
Stokes pulses is explored. 

CHOICE OF PULSE SHAPING 

The spectral phase of both pump and Stokes excitation 
fields is shaped simultaneously in such a way as to (1) 
suppress the strong-field induced level shifts, or (2) ini- 
tiate an AP process. In both cases, we aim at increasing 
the strength of the CARS signal. The important differ- 
ence between the two approaches lies in the two-photon 
pump-Stokes field. In scenario one, the spectral power 
densities at the desired Raman transition frequencies are 
maintained while the peak intensities of the pulses are 
reduced. Thus, the detrimental strong-field Stark shifts 
are minimized. In the second scenario, not only the field 
amplitudes of undesired frequencies but also those of the 
target Raman transitions are lower than in the TL case. 



Linear frequency chirping 

An equal frequency chirp applied to both pump and 
Stokes pulses results in the narrowing of the two-photon 
excitation spectrum around a constant frequency which 
is tuned to the frequency of the target vibrational tran- 
sition Qq. The spectral power density at that Raman 
transition is kept constant, whereas the peak power of 
the individual pulses can be substantially lowered. 

Linear chirping of the instantaneous frequency of pump 
and Stokes fields u!p s{t) — Wo(p,S) + is achieved by 
applying the parabolic spectral phase masks by means of 
two pulse shapers: 

V'p.slw) = - a;o(p,S))^- (1) 

The magnitude of the spectral chirp a' is related to the 
temporal chirp a — 1/a' [3T] . wop and wqs are the central 
frequencies of pump and Stokes pulses, respectively. An 
example of the temporal amplitude and phase before and 
after linear chirping is shown in Fig. [T|a). In comparison 
to a transform-limited pulse, the peak amplitude of the 
shaped pulse drops as the pulse is stretched in time in 
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FIG. 1: (Color online) Temporal profile of TL (dashed lines) 
and shaped (solid lines) femtosecond pulses with the elec- 
tric field amplitude (black single lines) and phase (red double 
lines), (a) Linearly chirped pulse with a' = 50,000 fs^, (b) 
sinusoidal phase modulation with A = 1.65, T = 334 fs. 

linear proportion with a'. The instantaneous frequency 
of the two-photon excitation spectrum of two linearly 
chirped pump and Stokes pulses with an equal chirp rate 
a remains constant, n{t) :— u!p{t) — to sit) — uJop — wqs, 
as shown in Fig. |2ja). Together with the two-photon 
field spectrum, we plot all possible Raman transitions 
in room-temperature iodine vapor within the accessible 
spectral bandwidth. Initial thermal population distribu- 
tion among the vibrational levels v = 0,1,2,... is taken 
into account. Vibrational resonances are broadened due 
to the thermal molecular rotation. We do not resolve the 
rotational spectrum in our experiments. The excitation 
line width is inversely proportional to the applied chirp 
a', whereas its frequency can be tuned by a variable time 
delay between pump and Stokes pulses. In CARS mi- 
croscopy this method is known as spectral focusing [llj . 
In our experiments, a' = 50, 000 fs^ is chosen to match 
the excitation line width with the rotational broadening 
of the vibrational transitions. 



Sinusoidal phase modulation 

Narrowing of the two-photon spectrum around the tar- 
get Raman frequency requires the spectral phase of 
pump and Stokes pulses to be identical around wop and 
ujQS = oJop — flo, respectively, as in the case of an equal 
frequency chirp described above. If the goal is to ex- 
cite a vibrational wave packet consisting of several eigen- 
states, the applied phase mask must be periodic in fre- 
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FIG. 2: (Color online) Raman spectrum at ujp — cos for TL 
pulses (black dashed line) and shaped pulses (red solid line). 
All Raman transitions (rotationally broadened into bands) 
and their strengths are indicated according to thermal Boltz- 
mann distribution at 100 °C: initial states v=0 (white area), 
v=l (grey area) and v—2 (black area), (a) Linearly chirped 
pulses with a' — 50,000 fs^, (b) sinusoidal phase modulation 
with A = 1.65, T = 334 fs. 

quency, with a period matching the vibrational period of 
a molecule [Hill [35]. 

Sinusoidal phase modulation represents one of the most 
popular methods of periodic shaping [32l [33l [36l [37] . In 

this case, 



(2) 



where T describes the modulation period and A its am- 
plitude, ujp and uis are the central modulation frequen- 
cies of pump and Stokes pulses, respectively. In the 
time domain, this spectral shaping results in a train of 
pulses, with each pulse being an exact replica of the ini- 
tial transform-limited pulse, and the time delay between 
the pulses defined by T (Fig. [IJb)). By matching or mis- 
matching the train period with the period of molecular vi- 
bration, the latter can be either enhanced or suppressed. 
Weak-field vibrational control by means of periodic pulse 
trains has been successfully demonstrated [BJ [21 [351 13E] ■ 
In the strong-field regime considered in this work, break- 
ing the initial pulse into a train of weaker pulses serves 
the purpose of suppressing the detrimental multi-photon 
processes prohibiting the transfer of molecules to the tar- 
get state. 

If an identical sinusoidal phase modulation is applied 
to both pump and Stokes excitation pulses, the two- 
photon spectrum shows periodic peaks with flat phase 
across them, see Fig.[2];b). Unlike the case of a frequency 
chirp, all three vibrational bands are now excited with 



the maximum possible amplitude (set by the available 
bandwidth) if the period T is a multiple integer of the 
vibrational period. The modulation amplitude A deter- 
mines the number of pulses in the pulse train, as well as 
the line width of the spectral peaks in the two-photon 
spectrum |3B]. Similarly to the parabolic phase shaping 
described above, we tune the value of A for the best cov- 
erage of the rotationally broadened Raman transitions 
in iodine. The central modulation frequencies Wp^s in 
Eq. [2] control the absolute position of the peaks in the 
two-photon spectrum. 



Pulse shaping for Adiabatic Passage 

In RAP, the instantaneous frequency of the excitation 
field is swept across the target transition frequency [3ll 
l39] . In the case of a Raman process, this corresponds 
to chirping the two-photon frequency. Applying linear 
frequency chirps with opposite signs to pump and Stokes 
fields ap = —0:5 = a, results in a linearly chirped two- 
photon field with an instantaneous frequency of Cl(t) — 
ujp(t) - ujsit) = ujQp - ujQS + 2at. 

The technique of STIRAP relies solely on a relative 
time delay between pump and Stokes pulses. Hence, no 
pulse shaping is necessary. The pulse trains needed for 
Piecewise STIRAP can be implemented via amplitude 
and phase shaping according to the method described in 
Ref. [27] (see also Fig.[5|. 



NUMERICAL ANALYSIS 

In our experimental detection scheme based on CARS, 
the signal from a single final state \vf) is proportional 
to the square of the coherence between that state and 
the initial vibrational state \vo), i.e. the square ampli- 
tude of |?;o)('i'/| in the density matrix. Since the rate 
of both collisional and rotational decoherence is negli- 
gible on the experimentally realized time scale of a few 
picoseconds, and the target population is small, the sig- 
nal magnitude is proportional to the population of the 
final state. We therefore base our numerical analysis on 
calculating the transfer of population from the ground 
vibrational state, predominantly populated at our exper- 
imental conditions, to higher vibrational states, driven 
by strong pump and Stokes excitation fields. For a co- 
herent wave packet consisting of several vibrational final 
states, the observed CARS signal oscillates as a function 
of the arrival time of probe pulses (see Results section 
below). Rather than model this time behavior, we as- 
sume that the time-averaged signal is proportional to the 
sum of the calculated populations. To calculate the re- 
distribution of the vibrational population as a result of 
the strong-field interaction, we solve the time-dependent 
Schrodinger equation in the eigenstate basis. The propa- 
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gation routine does not employ the rotating wave approx- 
immation and thus automatically accounts for AC Stark 
shifts. We take into account transitions between two elec- 
tronic manifolds of J2, X^'E'^^ and B^Hq^^- Both poten- 
tials are modeled as Morse oscillators [ID]. The elec- 
tronic transition dipole moment is set to an approximate 
value of 1 Debye, while the Franck-Condon factors for 
the B ^ X transitions are calculated numerically. While 
such modelling is not fully accurate for high vibrational 
states of iodine [H], it can provide a good qualitative 
description of the quantum dynamics. Quantitatively, in 
our simulations the strong-field effects arise at intensities 
somewhat lower than those observed experimentally. 

In the main set of calculations, we explore the sim- 
plest case of a strongly driven off-resonance Raman ex- 
citation. We avoid one-photon resonances, and set the 
central wavelength of the pump pulse to Ip = 800 nm 
whereas the Stokes pulse is varied between Is = 940 nm 
and 990 nm. Importantly, neither pump nor Stokes pho- 
tons have enough energy to drive a resonant transition 
between either \vq) or \vf) and the vibrational states of 
the B manifold. Various multi-photon interaction path- 
ways that couple \vo) to other electronic and vibrational 
states become intertwined in the strong-field limit. The 
population transfer is influenced by the time-dependent 
Stark shifts which move various vibrational levels in and 
out of resonance with different components of the Raman 
excitation spectrum. Nonlinear strong-field processes, 
most notably dynamic Stark shifts, make the perturba- 
tive approach inapplicable and lead to the observed sup- 
pression of CARS signal. 

We first focus on the pulse shapes designed to re- 
duce the peak intensity of the excitation pulses, while 
enabling the excitation of either a single vibrational level 
(quadratic spectral phase mask) or a superposition of lev- 
els (sinusoidal spectral phase mask). In the former case, 
a linear chirp with a' = 46, 000 fs^ stretches pump and 
Stokes pulses to ^ 1 ps duration. In the second case, the 
sinusoidal mask of Eq.([2]) with A = 1.23 and T = 312 fs 
leads to a pulse train consisting of five pulses separated 
in time by twice the oscillation period of a wave packet 
composed of the vibrational states adjacent to Vf = 10. 
The factor of 2 is introduced for an easier distinction of 
individual pulses in the resulting pulse train. 

Figure |3] demonstrates the dependence of the Raman 
transfer efficiency on laser intensity. Here, pump inten- 
sity {Ip) is fixed at 1.5 x 10^^ W/cm^ while Stokes inten- 
sity (Is) is scanned. The Stokes wavelength is set to 964 
nm corresponding to a strong two-photon resonance with 
Vf ^ 10. We compare the effect of 130-fs sin^ t/to-shapcd 
TL pump and Stokes pulses (full width at half maximum 
(FWHM) of the intensity profile) and that of the shaped 
pulses simulated with a numerical pulse shaper. As the 
experimental detection is carried out with a broadband 
probe and therefore does not offer frequency resolution, 
we plot the sum of populations in the target manifold. 
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FIG. 3; (Color online) Raman excitation probability as a 
function of Stokes intensity. The excitation schemes are: TL 
pulses (black solid line), pump and Stokes pulses are linearly 
chirped (blue dashed line), pump and Stokes pulses are shaped 
with a sinusoidal phase mask (red dotted line). 



V = 6. .13, rather than the populations of the individ- 
ual levels. This range is wide enough to include all the 
states covered by the excitation bandwidth in the pres- 
ence of strong-field Stark shifts. On the other hand, it 
does not include those states which can be populated in 
pump-pump and Stokes-Stokes Raman transitions in the 
presence of Stark shifts. 

The lower curve in Fig. [3] illustrates the main moti- 
vation of this work: the efficiency of the Raman exci- 
tation driven by high-intensity transform-limited pulses 
saturates, setting an upper bound on the magnitude of 
the CARS signal. The saturation stems from an os- 
cillatory behavior of the individual target populations, 
which takes over the linear (with Is) growth of popu- 
lation anticipated in the weak-field limit. We associate 
these population oscillations with an interplay of popula- 
tion transfer with multiple AC Stark shifts which dynam- 
ically change the instantaneous spectrum of the molecu- 
lar states dressed by the strong laser field. We note that 
on the scale of laser intensities considered in this work, 
the oscillation amplitudes are typically of the order of 1- 
2% of the total population, which explains the low total 
efficiency of the Raman excitation plotted in Fig.|3] If the 
two-photon field is resonant with a Raman transition to 
a certain vibrational state then the linear frequency 
chirping of pump and Stokes pulses provides the largest 
suppression of the strong-field saturation (dashed curve 
in Fig. |3|. Raman excitation by a train of pump/Stokes 
pulse pairs (sinusoidal spectral mask) showed an inter- 
mediate level of saturation (dotted curve in Fig. [3| . We 
attribute this difference in performance to the fact that 
for a given pulse energy, frequency chirping results in the 
lowest peak intensity and therefore better suppression of 
the detrimental strong-field effects discussed above. 

Further insight can be gained by calculating the trans- 
fer efficiency as a function of the central Stokes frequency, 
ujQS- The results are shown in Fig. |4j In the weak-field 
regime (left panel), the transition strength is defined by 
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FIG. 4: (Color online) Raman excitation probability as a 
function of Stokes wavelength. Left; both pump and Stokes 
intensities are set to 10^^ W/cm^. Right: both intensities are 
set to 10^"^ W/cm^. Line assignment is the same as in Fig.js] 

the corresponding resonant component of the two-photon 
pump-Stokes spectrum. Since the peak spectral ampli- 
tudes are equal for all three pulse shapes, the population 
transfer to a single state \vf) is equally efficient when 
ujQp — loqs = {Ef — Eoj/h, where Ej. is the energy of 
the fc-th vibrational state. In this figure, the Raman 
resonance with Vf — 10 lies at Is = 964 nm. Raman 
transitions to the neighboring = 9 and Vf = 11 are 
much weaker due to oscillations in the two-photon cou- 
phng strength. For the frequency chirped pulses, the 
two-photon spectrum is narrower than in the case of TL 
pulses (see Fig. [2|a)), resulting in a faster drop of the 
transfer efficiency away from the resonance. The two- 
photon spectrum of the pump-Stokes pulse train has a 
much wider maximum [44 . 

The right panel of Fig. |4] shows the calculated trans- 
fer efficiency in the strong-field limit. Dynamic Stark 
shifts modify the line shape and lead to a significantly 
lower population of the target vibrational manifold. The 
vibrational excitation by a train of pulses is suppressed 
to a much lower degree. Indeed, the intensity of pulses 
in the train is lower than that of the original femtosec- 
ond pulse and the strong-field effects are correspondingly 
weaker. Finally, chirped pulses correspond to the highest 
population transfer. The absence of a visible line broad- 
ening in this case points to the complete elimination of 
the detrimental strong-field effects. 

Our calculations do not indicate that the adiabatic 
techniques of two-photon RAP EH El], STIRAP [7 or 
piecewise STIRAP [32l [33] are feasible. Indeed, AP is 
usually used to improve the robustness of an already 
high transfer rate. The left panel of Fig. [5] investigates 
the strong-field transfer for the three adiabatic schemes. 
The intensity profiles are shown in the right panel with 
wavelengths and intensities similar to those in Fig. [3j 
For RAP, we show the transfer efficiency with pump 
and Stokes fields chirped in the opposite directions with 
\a'\ = 46,000 fs^. This shaping results in a frequency 
chirp of the two-photon pump-Stokes field. For STI- 
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FIG. 5: (Color online) Left: Raman excitation probability 
as a function of Stokes intensity. The excitation schemes 
are: RAP (purple dashed line), STIRAP (orange dotted line), 
piecewise STIRAP (green dash-dotted line) and TL pulses 
(black solid line). Right: Schematic view of the intensity pro- 
files of pump (blue dashed line) and Stokes (red solid line) 
pulses used in the calculations of RAP (top), STIRAP (mid- 
dle), and piecewise STIRAP (bottom). 

RAP, transform-limited Stokes pulses precede transform- 
limited pump pulses by 130 fs. For piecewise STIRAP, we 
consider two pulse trains obtained via sinusoidal spectral 
phase modulation as described above and shifted relative 
to each other in the counter-intuitive order by the pulse 
train period [32 ES]. In all three cases, the two-photon 
spectrum is lower than that of a TL pulse. Hence, the 
transfer efficiency at small intensities is lowered by the 
applied pulse shaping. At high intensities, the transfer 
efficiency is only slightly higher than that of transform- 
limited pulses, but still is significantly lower than that 
achieved by the non-adiabatic techniques in Fig. [3| 

Finally, we compare the results of the off-resonance 
excitation regime, with the two resonant cases, where a 
single-photon electronic excitation from either wq (case 
1, /p = 573 nm. Is = 800 nm) or (case 2, Ip = 645 
nm. Is = 800 nm) is energetically allowed. The choice 
of these wavelengths is based on our experimental set- 
tings and on the previous studies [2S1[21]. In both cases, 
one-photon coupling to the excited electronic state at the 
pump frequency significantly complicates the population 
dynamics. The dependance of the transfer efficiency on 
the wavelength is highly irregular due to the interplay of 
time-dependent AC Stark shifts and a single-photon es- 
cape of population into the B manifold. In the first case, 
the overall transfer efficiency at high intensities is an or- 
der of magnitude lower than that in the far off-resonance 
arrangement considered in the main set of calculations. 
In the second case (previously studied in Ref. [28 ] [29 ] ). 
the efficiency of the transfer can reach up to a few per- 
cent, i.e. slightly higher than that in the far off-resonance 
arrangement; however, the dependance on the laser pa- 
rameters remains erratic. Here, the efficiencies observed 
in our calculations are substantially lower than those re- 
ported in Refs. [2H1[5^. In our calculation, the target 
population transfer into w/ is hindered by one-photon 
coupling to the excited electronic state at the pump fre- 
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FIG. 6: (Color online) Experimental setup: (a) Colinear 
CARS geometry with spatial filtering, (b) Use of polarizing 
optics and spectral filtering to separate CARS signal from the 
excitation pulses. 



quency. 



EXPERIMENTAL STUDY 



Experimental setup 

A Ti: sapphire-based laser system (SpitFire Pro, 
Spectra-Physics) produces 2 mJ 130 fs pulses at 800 nm 
and 1 kHz repetition rate. These pulses serve as pump, 
whereas one optical parametric amplifier (TOPAS, Light 
Conversion) generates Stokes pulses at a wavelength of 
973 nm and another OPA (OPA-800C, Spectra Physics) 
generates probe pulses at 578 nm. 

All three pulses are aligned coUinearly and are spa- 
tially overlapped inside a vapor cell filled with I2 at vari- 
able temperature, see Fig. [6]ja). At 100°C, used in this 
work, the vapor pressure is 43 Torr [32] • Collinear ge- 
ometry is implemented to increase the interaction length 
and reduce the effect of spatial averaging over the Gaus- 
sian beam profiles. The latter is also achieved by mak- 
ing the focal size of the probe beam smaller than the 
size of the other beams. With a 30 cm focusing lens, 
peak intensities of up to 5 x 10^'^ W/cm^ for pump and 
Stokes pulses and 2 x 10^^ W/cm^ for probe pulses are 
obtained. Spatial filtering in the detection channel is 
used to eliminate a majority of the fluorescence back- 
ground. In order to filter the output anti-Stokes beam 
from the incident beams, we use polarization and spec- 
tral filtering, as shown in Fig. ^h) . Probe pulses are lin- 
early polarized in the orthogonal direction to pump and 
Stokes pulses. Hence, the anti-Stokes polarization does 
not coincide with that of the strong excitation fields, en- 



abling one to block the latter while passing through part 
of the CARS signal [43]. The filtered signal is sent to 
a spectrometer equipped with a cooled (— 40°C) charge- 
coupled device (CCD) camera. Two delay lines are used 
to vary the relative arrival time of pump. Stokes and 
probe pulses. Pump and Stokes pulses are shaped by 
two separate pulse shapers implemented in a standard 
4/-geometry [2]. 



Results 

To evaluate the effect of strong excitation fields, we 
start by detecting the output spectrum in the absence of 
probe pulses. The observed spectrum, shown in Fig.[7]^b), 
is very broad and covers wavelengths above 400 nm. It 
corresponds to spontaneous emission following the pro- 
cess of two-photon absorption (TPA) in which any com- 
bination of two photons from the high intensity pump 
and Stokes fields are absorbed. The 400 nm cut-off is 
determined by the maximum energy from the absorp- 
tion of two 800 nm pump photons. If only pump or 
Stokes photons are present, the fluorescence signal shows 
a quadratic dependency on the intensity of the corre- 
sponding laser beam. The fluorescence signal stemming 
from TPA of one pump and one Stokes photon shows a 
linear dependence on the intensity of both laser beams. 
It provides the main contribution to TPA (Fig. [7][a)). 
This is expected since the two-photon excitation ampli- 
tude scales as {Ep + EsY = Ep + E'g + 2EpEs, resulting 
in a four-fold enhancement of the intensity dependence 
on Ipis (Here Ep{Ip) and Es (Is) is the electric field 
amplitude (intensity) of pump and Stokes pulses, respec- 
tively). The exact proportionality factor is a function 
of the frequency dependent Frank-Condon factors. The 
strong peak at 678 nm is due to the non-resonant four- 
wave mixing process with two pump photons and one 
Stokes photon at frequency 2wp — ujs- 

Strong two-photon fluorescence, comparable in magni- 
tude to the detected CARS signal even after both spatial 
and polarization filters have been applied (Fig. [Tjja)), re- 
flects a high degree of two-photon coupling responsible 
for the detrimental Stark shifts. Suppressing the TPA- 
induced fluorescence with pulse shaping can be consid- 
ered as an indirect evidence of lowering the influence of 
strong-field effects on the target Raman excitation. Both, 
the linear frequency chirp and the sinusoidal phase mod- 
ulation, drastically reduce the fluorescence background 
due to the narrowing of the two-photon spectrum dis- 
cussed earlier in the text. The effect is shown in Fig. |8| 
As expected, linear chirping achieves better suppression 
of fluorescence, since the two-photon spectrum has fewer 
spectral components. In contrast to the suppressed TPA, 
the intensity of the CARS signal is rising, confirming the 
above argument about the mechanism of the strong-field 
induced saturation. 
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FIG. 7: (Color online) Fluorescence spectrum of iodine due 
to strong pump and Stokes excitation, (a) Detected signal in 
the wavelength range from 482 nm to 522 nm stemming from 
TPA of Stokes (blue dashed), TPA of pump (green dotted) 
and TPA of mixed pump and Stokes (black solid). CARS 
signal is shown on top of the fluorescence curve (red shaded 
area) . (b) Fluorescence spectrum in the range from 350nm to 
700nm. 
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FIG. 8: (Color online) Detected spectrum with strong pump 
and Stokes fields and a weak probe field (red dotted line), 
output spectrum in the absence of probe pulses (blue dashed 
line), CARS signal derived as the difference of the two spec- 
tra (black solid line); Pump and Stokes pulses are both (a) 
transform-limited, (b) modulated with a sinusoidal phase, (c) 
linearly chirped. 

In Fig. [9] we plot the detected CARS signal obtained 
by subtracting the fluorescence background from the 
measured spectrum. We note that in this spectral re- 
gion (i.e. around ujp — ujs + ujpr), the signal is free of 
non-resonant background since probe pulses are delayed 
by about 25 ps with respect to the temporal overlapping 
pump and Stokes pulses. The maximum available energy 
of the excitation pulses corresponds to peak intensities of 
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FIG. 9: (Color online) Comparison of CARS intensity with 
different excitation schemes: Pump and Stokes pulses are 
both transform-limited (black solid line), linearly chirped 
(blue dashed line) or shaped with a sinusoidal phase mask 
(red dotted line), (a) CARS spectrum, (b) CARS signal as 
a function of the probe time delay and (c) the corresponding 
Fourier spectra. 



ip,S = 5x 10^'^ W/cm^ of the unshaped TL pulses. As can 
be seen in plots (a) and (b) of Fig. [Oj both phase masks 
(quadratic and sinusoidal) result in an almost identical 
signal enhancement by about 50% in comparison to the 
case of the unshaped excitation. 

The dependence of the CARS signal on the probe ar- 
rival time is plotted in Fig. [9f^b) . Aside form the similar 
difference in the absolute signal strength, we note that 
the oscillations of the excited vibrational wave packet, 
clearly seen in the case of the TL excitation, are preserved 
in the case of the sinusoidal phase modulation and sup- 
pressed in the case of the linear frequency chirping. The 
Fourier spectra in Fig. [ojjc) reveal a peak at 200 cm~^. 
This corresponds to a vibrational wave packet consisting 
of Vf — 10, 11 and 12 with an oscillation period of 167 fs. 
The peak is suppressed when frequency chirped pulses are 
used. The effect is different from a well known feature of 
the weak-field regime, where the excitation amplitude is 
linearly proportional to the amplitude of the two-photon 
field at resonant frequencies. Although the sinusoidal 
phase modulation does not change the spectral amplitude 
at the frequencies of the available Raman transitions (see 
Fig. [2];b)), it resuhs in a substantial signal gain while 
allowing us to excite a coherent wave packet. On the 
other hand, the observed signal enhancement with the 
frequency chirped pump and Stokes pulses is accompa- 
nied by the two-photon spectrum narrowing (Fig. [2][a)), 
and a consequent loss of the wave packet oscillations. 

By scanning the intensity of the Stokes beam in Fig.fTO} 
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FIG. 10: (Color online) CARS intensity as a function of 
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we observe the effect of pulse shaping above approxi- 
mately 1.5 X 10^'^W/cm^ when the signal from the un- 
shaped transform-limited excitation starts to slow down. 
We attribute this onset of saturation to the deleteri- 
ous strong-field effects, which are suppressed when the 
pulses are shaped. In the latter case, the signal keeps 
growing almost linearly exceeding the unshaped limit by 
about 50%. Our numerical analysis suggests that fre- 
quency chirping should provide the highest CARS signal 
(see Fig. |3]), whereas in the experiment, the two spectral 
masks result in a similar performance. In the simula- 
tions, we set the Stokes wavelength to match the peak 
of the two-photon spectrum with a Raman resonance. 
In the experiment, however, the wavelengths cannot be 
changed easily and the frequencies are matched by in- 
troducing a delay between pump and Stokes pulses (see 
also next paragraph) . This delay reduces the two-photon 
spectral amplitude, since it is shifted from the maximum 
of the two-photon spectrum, thus lowering the signal. 
In the case of the sinusoidal phase modulation, the ex- 
act overlap of the two-photon spectrum with the Raman 
transition frequencies is easily maintained by adjusting 
the central modulation frequencies ojp^s of Eq. [2] with 
pump and Stokes pulse shapers. 

Following our numerical analysis, we examine the tech- 
niques of RAP and STIRAP in both the off-resonance 
and near-resonance cases. In the off-resonant case, i.e. 
with the same choice of pump and Stokes wavelengths as 
before { Xp = 800 nm, A5 = 973 nm), we find no evidence 
of the AP enhancement of CARS signal. As required for 
rapid adiabatic passage, we chirp pump and Stokes fields 
in opposite directions with \a'\ = 50, 000 fs^, and observe 
CARS intensity always decreasing below the TL limit. 
This agrees with the perturbative, rather than adiabatic, 
description of the Raman process. As discussed earlier 
in the text, opposite-sign chirping of pump and Stokes 
pulses results in a lower two-photon intensity which is 
not sufficient for initiating AP. 



In order to test the non-perturbative regime necessary 
for AP, we shift the central excitation wavelengths to 
Ip — 645 nm and Is = 800 nm. This choice of wave- 
lengths rephcates that of Refs. [21 [gS]. In Fig. [TTj^a), 
the spectra of the excitation pulses are chirped linearly 
in opposite directions (RAP) with two different chirp am- 
plitudes a'. CARS intensity is plotted as a function of 
the time delay between the pulses. We find no signatures 
of RAP (i.e. signal increase with frequency chirping at 
zero time delay) or STIRAP (i.e. signal increase with 
time delay at zero frequency chirp), and always record 
the strongest CARS signal with transform-limited pump 
and Stokes pulses perfectly overlapping in time. From 
this observation, we conclude that no adiabatic evolution 
takes place. 

In Fig. 11 [h,c), we apply the successful same-sign fre- 
quency chirping, found in the off-resonance case, to the 
near-resonance interaction scheme. As discussed earlier, 
the two-photon spectrum becomes narrower and, depend- 
ing on the exact pump-Stokes time delay, only one Ra- 
man transition is covered at a time (c.f. Fig. [2]ja)). This 
is reflected by the peaks in the CARS intensity at cer- 
tain time delays seen in Fig.[lT];b). Fig.jnjc) shows that 
for all beam intensities achievable in the experiment, the 
overall CARS intensity is lower than that for the un- 
shaped pulses. Unlike the off-resonance case, pulse shap- 
ing does not improve the overall CARS efficiency. We 
attribute the difference to complicated dynamics due to 
one-photon coupling to the excited electronic state. 



SUMMARY AND DISCUSSION 

In this experimental and theoretical work, we inves- 
tigated the utility of femtosecond pulse shaping to in- 
crease the efficiency of CARS in the strong-field excita- 
tion regime. We showed that when molecular eigenstates 
are significantly modified by the applied laser fields, the 
efficiency of exciting molecular vibration is no longer 
achieved with unshaped transform-limited pulses. We 
analyzed two qualitatively different approaches to en- 
hancing the magnitude of non-resonant vibrational ex- 
citation with shaped pulses. 

In the first approach, pulse shaping was used to sup- 
press strong-field effects, such as AC Stark shifts, while 
preserving the amplitude of the two-photon field at the 
frequency of Raman resonances. Numerically, we showed 
that linear frequency chirping of strong pump and Stokes 
pulses can provide several times better excitation effi- 
ciency than that achieved with transform-limited pulses. 
The efficiency of Raman excitation with sinusoidal phase 
modulation showed an intermediate result. In the ex- 
periment, both techniques improved the strength of the 
observed CARS signal by 50% with respect to the TL 
case. Sinusoidal phase modulation resulted in a substan- 
tial signal gain while allowing us to excite a coherent 
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efficiency of Raman excitation. We associate this lack of 
success with two factors. Our numerical analysis showed 
that the complexity of the molecular spectrum prevents 
reaching the AP regime which is typically manifested by 
a high amount of population transferred to the target 
state or set of states. From the experimental point of 
view, adiabatic passage requires shaping the excitation 
pulses in a way which reduces the two-photon spectral 
power density, making the available laser power insuffi- 
cient for establishing an adiabatic evolution. 

The authors would like to thank Vladimir Malinovsky 
and Sergey Zhdanovich for valuable discussion. This 
work has been supported by the DTRA grant HDTRAl- 
09-1-0021. 
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FIG. 11: (Color online) Near-resonance CARS, (a) CARS 
intensity as a function of pump-Stokes time delay for TL 
pulses (black solid line) and oppositely chirped pulses with 
a's^p = ±10,000 fs^ (blue dashed line) and ag ^ = ±25,000 
fs^ (red dotted line), (b) CARS intensity as a function of 
pump-Stokes time delay and (c) as a function of pump inten- 
sity with Is = Ax lO^^W/cm^ for TL pulses (black solid line) 
and same-sign linearly chirped pulses with a'g p = —10, 000 fs^ 
(blue dashed line) and a's^p = —25, 000 fs^ (red dotted line). 

vibrational wave packet, which may prove important in 
strong-field spectroscopic applications. 

Multiple sets of pump and Stokes wavelengths were 
tested. Our main study was devoted to the far off- 
resonance case, where the dependence of vibrational ex- 
citation on the field parameters allows reasonably simple 
interpretation. In the near-resonance case, when either 
the ground or the target state were coupled to the ex- 
cited electronic state, the population dynamics were er- 
ratic. In spite of the higher two-photon Raman matrix 
elements, the excitation efficiency in all studied near- 
resonance cases was either negligibly higher or even sub- 
stantially lower (depending on the choice of wavelengths) 
than in the far off-resonance case. 

In the second approach, we put to test three adia- 
batic passage scenarios which excel in a robust and ef- 
ficient population transfer in simple two- or three-level 
atomic systems - RAP, STIRAP and piecewise STIRAP 
[HIISII^. All AP schemes studied in this work failed to 
establish an adiabatic evolution and did not improve the 
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